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Smoothing is an estimation method whereby a classical state (probability distribution for classical 
variables) at a given time is conditioned on all-time (both earlier and later) observations. Here we 
define a smoothed quantum state for a partially monitored open quantum system, conditioned on 
an all-time monitoring-derived record. We calculate the smoothed distribution for a hypothetical 
unobserved record which, when added to the real record, would complete the monitoring, yielding 
a pure-state “quantum trajectory”. Averaging the pure state over this smoothed distribution yields 
the (mixed) smoothed quantum state. We study how the choice of actual unravelling affects the 
purity increase over that of the conventional (hltered) state conditioned only on the past record. 

PACS numbers: 03.65.Ta, 03.65.Yz, 42.50.Dv, 42.50.Lc 


Estimation theory is used to assign values to param¬ 
eters of interest, whose true values are unknown, using 
the available data. These parameters may evolve dynam¬ 
ically, and new data may arrive dynamically, through a 
continuous measurement process. Estimation in this in¬ 
stance is non trivial because there can be noise associ¬ 
ated with the measurement, noise affecting directly the 
dynamical system due to its environment, and initial un¬ 
certainty in the parameters. Optimal estimation theory 
can be formulated using the Bayesian approach to statis¬ 
tics, whereby the observer’s knowledge of the parameters 
is described by a conditional probability distribution pc , 
also known as a Bayesian state. If this state is conditional 
on measurements at earlier times, it is called a filtered 
state pF, while if it is conditional on all-time (both ear¬ 
lier and later) measurements, it is called a smoothed state 
ps- Smoothing uses more complete information than fil¬ 
tering, and so typically delivers a probability distribution 
that is purer (that is, having less entropy). 

In the flourishing and important area of quantum esti¬ 
mation theory HHH], much has been learnt from classical 
estimation theory. The analogy between quantum states 
and classical Bayesian states has been fruitful even in 
quantum foundations [uns]. In particular, the stochas¬ 
tically evolving conditioned state pc of an open quantum 
state, as introduced by physicists USM] and applied in 
quantum control p1 [22ll29] . is now understood to be anal¬ 
ogous to the classical filtered state ISHEDHSl pF (and 
so for clarity we write it as pf)- However, the situation 
is very different regarding smoothing. 

The term “quantum smoothing” was introduced by 
Tsang [35l [36] in 2009 to mean smoothed estimation of 
classical parameters that affect the evolution of a quan¬ 
tum system, using the results of measurements on that 
system. It has been shown to be useful to the prob¬ 
lem of estimating a stochastically varying optical phase 
using the all-time photocurrent record, both theoreti¬ 
cally [33138] and experimentally [39l|40]. It has also been 
applied to the problem of estimating an unknown result 


from a measurement on a quantum system at one time, 
using records obtained both before and after that time, 
again both theoretically [41] and experimentally [42] . 
However, none of the above define a quantum smoothed 
state—that is, a positive operator ps that is analogous 
to a Bayesian smoothed state ps- There is a good reason 
for this lack, which is most easily stated in the Heisen¬ 
berg picture [35]: quantum operators for a system at a 
given time commute with operators representing the re¬ 
sults of earlier measurements on that system, but do not 
commute with operators representing the results of later 
measurements on that system [T] . 

In this Letter we show that there is a situation in which 
it is possible to define quantum state smoothing, produc¬ 
ing a positive state ps conditioned on both earlier and 
later results. The situation is that of open quantum 
systems under partial observation, which is the typical 
situation experimentally [26[ [23 [29]. The system has 
couplings to several baths (all assumed Markovian). An 
experimenter, Alice, is able to monitor some of them, 
yielding the record O she observes. Other baths are not 
monitored by her, but hypothetically they could be mon¬ 
itored by another party, yielding results U unobserved 
by Alice. The ‘true’ state px conditioned on both ob¬ 
served (O) and unobserved (U) records would be pure, 
while that conditioned only on O is mixed. The crucial 
point is that the record U, comprising classical variables 
(c-numbers), can be estimated by applying smoothing to 
the record O, and in this way Alice can obtain a smoothed 
quantum state ps- As in the classical case, this is typi¬ 
cally purer than pF, and a better approximation to px- 

We first review the necessary theoretical background. 
Then we explain how our smoothed quantum state ps 
is quite different from the “past quantum state” (actu¬ 
ally a pair of operators) introduced in Ref. [41]. We also 
show that our approach subsumes classical state estima¬ 
tion by smoothing, also known as the hidden Markov 
model (HMM) technique [43]. A HMM is applicable to 
quantum systems only when they are effectively classical 
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(i.e. always diagonal in a particular basis) as in Refs. [44]- 
[46] , We apply our method to a genuinely quantum sys¬ 
tem (i.e. one which is not diagonal in a fixed basis)—a co¬ 
herently driven two-level atom, the radiation from which 
is partly observed. We take the known record to be gen¬ 
erated by homodyne detection, and the unknown record 
to be that corresponding to photon absorption, as this 
is the most intuitive picture of what happens to photons 
that are lost into the laboratory surroundings. These 
lost photons result in impurity in the standard (filtered) 
conditioned system state. We show that our smoothing 
technique can, on average, eliminate up to 26% of this im¬ 
purity. Our investigations shed light on how well we can 
know the trajectory of a partially observed open quan¬ 
tum systems, and the relation between the quantum and 
classical versions of state smoothing. 

Types of Estimation .— Consider a classical dynamical 
system described by parameters (bold font indicates a 
vector of parameters) which is monitored to yield a noisy 
output at each time, r^. We denote a measurement reeord 
Hq = {rt : t G where Q C [to,T] is typically some 
finite time interval. Bayesian estimation involves data 
processing to infer the conditional classical state 

pRsj(xr) = Pr[x*™® = x^|Rn; po], (1) 


where po describes the a priori statistics of x at the initial 
time to- It is also useful to define the unnormalized state 


PRs2(Xr) 


PRsi(Xr) 


p(R-n|po) 

Post(I^nlpo) 


(X p(Rn,x^|po). (2) 


Here the p(Rq|po) is the actual distribution for while 
post(Ro|po) is an ‘ostensible’ distribution for it—it is 
positive and normalized, but is otherwise arbitrary and 
does not depend on Xt m- 

There are three types of estimation worth distinguish¬ 
ing [48l [49]: filtering, retro-filtering (as we call it), 
and smoothing (see Fig. [^. If—as in feedback control 
problems—for the time of interest r there is only access 
to earlier results, the optimal protocol is 

filtering: pf(xt-) = p^ (^r). If there is access only to 

later results, R^- = R[r,T): fhe optimal protocol is retro- 
filtering: pr(xt-) = pg (xr)- As its name implies, this is 
simply the time-reverse to filtering, but starting with an 
uninformative final state p(xf) oc 1. Finally, if the all- 
time record ^ is available, with to < r < T, 

then all the information can be utilised, by the technique 
of smoothing: Ps(xt-) = p^(xt-). This combines filtering 
and retrofiltering, using unnormalized states [36] : 


/ X ^ Pr(x^)pf(x^) 

/ dx; PR (x; ) pf (x; ) ■ 


(3) 


Here one of the states (most conveniently the retrofiltered 
one) is defined using an uninformative prior, po oc 1, to 
prevent double counting of the a priori information. 



FIG. 1: Classical estimation classes depending on the mea¬ 
surement record considered relative to r, time at which the 
signal is to be estimated. (Adapted from [35]). 


Quantum analogues of (retro)filtering. — An extension 
of these results to quantum mechanics has been done par¬ 
tially. Quantum trajectory theory m is the analogue of 
classical state filtering. A quantum trajectory describes 
the path (in “density operator space”) of the state of the 
quantum system through time, conditioned on the mea¬ 
surement result Tf in each infinitesimal interval [t,t-\-dt). 
Note that the path may be continuous (quantum diffu¬ 
sion) or discontinuous (quantum jumps) [T] [16]. This 
process is described by a set (indexed by r^) of measure¬ 
ment operations (completely positive maps) A4r^, that 
evolve the state forward in time ST]: 

Pf {t + dt) = MrtpF (t) . (4) 

Starting with p{to) = po^ this procedure generates the 
state conditioned on the whole past record: Pf(t) = 
p^ (t). This is an unnormalized state (as indicated by 
the tilde), analogous to Eq. (H- That is, the normalized 
version Pf(t) generates the correct filtered probability 
distribution pF(^r) for any system observable X^-, while 

Tr[pF(r)]post(kr|/)o) = p(kr|po)- (5) 

The corresponding analogue for Bayesian state retro- 
filtering has been set out in [35]; it is the solution of the 
adjoint of equation 0. 

Enlt) = ER(t + dt). (6) 

In this case the effeet operator evolves backwards from 
the final uninformative effect E(T) = I towards Eyi^t) = 
E^ (r), conditioned on the record R^- in the future of 
r. This solution E^ (r) determines the statistics of R^-: 

Tr[-^R('r)PT]post(ST|PT) = p(StIPt)- (7) 

Quantum smoothing? — A naive approach to construct 
a quantum smoothed state, given the quantum analogues 
of filtering pF and retrofiltering Er, would be to com¬ 
bine them directly as in equation Eq. 0 so that Ps{t) oc 
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FIG. 2: The quantum state smoothing problem: to best ap¬ 
proximate the unknown true state of a quantum system, con¬ 
ditioned on both observed ( O) and unobserved ( U) records, 
given access only to O. This requires one to estimate tjt up 
to time t using the full record for O (before and after t). 


PF('r)^R(r). However, as pointed out in Ref. [35], the re¬ 
sult is not in general Hermitian or (even if symmetrized) 
positive semi-definite. Therefore it cannot correspond to 
a physical state. As discussed in the introduction, there 
is a deep reason for this, which is why Tsang gave quan¬ 
tum smoothing the restricted meaning of estimating an 
external classical parameter x. 

The filtered state Pf{^) and the retrofiltered effect 
are sufficient to best estimate, from the all-time 
record the result of any measurement performed on 
the system at time r. For this reason, Gammelmark et 
al. declared the pair E = (pp, ^r) to be “the past quan¬ 
tum state” m- They did not define a quantum state in 
the usual sense {i.e. a density operator) combining the 
future and past information. While our notion of quan¬ 
tum state smoothing also makes use of filtered states and 
retrofiltered effects, it is quite different in that: it applies 
only to partially observed open quantum systems; it de¬ 
fines a quantum state in the usual sense; and it can be 
compared directly to quantum state filtering by measures 
such as purity, similarly to the classical case. 

Quantum state smoothing .— The key idea for quantum 
state smoothing is illustrated in Fig.[^ Consider an open 
quantum system with two groups of output channels (5, 
c). An observer Alice monitors the first group, 6, yield¬ 
ing the all-time measurement record O . A hypothetical 
observer monitors the second group c, yielding a record 
t? unobserved by Alice. The ‘true’ state of the system 
pT{t) = conditioned on both measurement 

records. If po is pure then Pt(^) will be pure for all times; 
no extra conditioning could possibly give a better (more 
pure) state. However Alice’s conditional state, calculated 
in the conventional way (filtering). 


will be mixed, because of the averaging over the unob¬ 
served record—here E^|^[X] means the expected value 
of X, averaged over A, for a given B. Note that this 
averaging does not have to be done explicitly—it is im¬ 
plicit in the quantum trajectory theory as in Eq. and 
is independent of how the channels c are monitored. 

The crucial insight is that Alice can do better, when 
averaging over tj^, by using information in the future of 
t, to define a positive-definite smoothed quantum state 

psit) = (9) 

tit 


Here ps(E) = p^(E) = = E|b,po] is 

the probability distribution for the unobserved record 
prior to t, obtained by smoothing from the all-time ob¬ 
served record Note that ps{^t) yields exactly the 
same type of information as the “past quantum state” 
of Ref. [41], except that it is more general — it speci¬ 
fies the probability of a continuous monitoring record tj, 
not just a result of a measurement at one point in time, 
and is conditioned on another record, O , covering the 
same time interval, not just records strictly earlier and 
strictly later. Unlike the classical stochastic process Xt 
considered previously by Tsang [35], the record tjt is of 
quantum origin—its statistics are undefined without a 
quantum system. Nevertheless, it is still a time-series of 
c-numbers with well-defined statistics and so there is no 
conceptual problem in applying his theory of quantum 
smoothing to obtain ps(Ut), and thereby ps{t). 

We now show how to calculate Eq. Eor definiteness 
and simplicity, we consider a single channel (6) yielding 
homodyne photocurrent yt and a single channel (c) yield¬ 
ing an unobserved photon count These processes are 
related to the dynamics of the quantum system via the 
joint measurement operation Mnt,yt defined such that 

■^Vt = I]nt=o PoBt{nt\Vt)Mnt,vt, for a convenient choice 

of post{nt\yt) mi- By standard techniques [T], M.rit,yt 
lets us generate a typical sample of the all-time records 
0 = 1 and ^ For all but one purpose 

(see below), the latter is irrelevant, but using the former 
we calculate the filtered state (t) from Eq. Q and the 

retrofiltered effect (t) from Eq. , with A4y^ in place 


of M-vt [S3]. Next, we generate a large ensemble (^ost of 
random samples of ^ = ^ ^ according to the ostensible 
distribution post(^tl^)- For each sample we calculate 
an associated pure state conditioned on both 

records, from Eq. Q with Mm^yt place of Afr* [54] . 
Elementary manipulation of probabilities m gives 

Ps{%) = pAt|T) « Using the 

equations for multiple channels corresponding to Eq. 0> 


( 8 ) p{t\%M) ( 10 ) 


PF(f) = P^p) - 
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and to Eq. (§, 

Post(^l^) = P(^l^), (11) 


we 


finally obtain, from Eq. ©, 


Ps{t) oc y^post(^li^) X 

( 12 ) 

We can approximate this weighted average over all pos¬ 
sible unobserved records using the ensemble (^ost drawn 
from the appropriate ostensible distribution, as discussed 
in the preceding paragraph. This is the method we use 
below to find the smoothed quantum state. 

Example .— Consider a two-level atom, with driving 
Hamiltonian H = {Q/2)ax in the interaction frame, 
and radiative damping described by a Lindblad opera¬ 
tor y/qd- [ 1 ]. We take a fraction rj of the fluorescence 
to be observed by homodyne detection, so b = 

The remainder is absorbed by the environment, which 
we model as an unobserved record of photon counts, as 
discussed above, with c = y^ 7(1 — 77 )d_. Eor a fixed 
we can compare ps with pp on the interval [ 0 , T]. At the 
final time ps{T) = Pf(^) because there is no more fu¬ 
ture record Ip to give extra information to ps{T). Also, 
we take the initial state to be pure, po = | 1 ) (ll^ which 
guarantees that px is pure and that ps(0) = Pf(0). 

To evaluate the advantage gained by smoothing over 
Altering, we use the purity. 


p[pcm=npUt)], 


(13) 


where pc could be either pp or ps- If (as is the case in 
simulations) we know the ‘true’ unobserved record 
we can also calculate the fidelity of the conditioned state 
to the true state pT(t) = 

F[pT{t),pc{t)] = Tr[pT(t)pc(t)]. (14) 

It is easy to show m that these measures are related by 
E {Pipcim = E {F[pT{t), pc(i)]} (15) 


where the ensemble averages here are over the actual dis¬ 
tributions for and r . 

In Eig. I^a) we show typical trajectories, for Y- 
homodyne [4> = 7r/2] for a randomly generated true state 
Pt(^) = featuring one jump at t 1.8. In this, 

case, Px, Ps, and pp are all confined to the Y-Z plane of 
the Bloch sphere, as shown. We plot Eqs. ([T^-([T^ in 
Eig. [^b) and (c) respectively. It is notable from (b) that 
Ps anticipates the jump in px and its uncertainty about 
the timing of the jump leads to a lower purity in the 
region of the jump than the non-anticipating pp. Sim¬ 
ilarly, (c) shows that the fldelity of ps to px decreases 
below that of pp prior to the jump, but is higher after 
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FIG. 3: [Color online.] (a) Trajectories in the Bloch sphere 
for our model system with Q = 20, 7 = 1, and p = 10/11. The 
states shown are pp (filtered, blue), ps (smoothed, purple) and 
Pt (‘true’, green) for a case where the ‘true’ record includes a 
jump. We also plot the purities (b) and fidelities with px (c) 
of these pp and ps- The purities for a record with no jump are 
shown in (d). To compute Pg we average over an ensemble of 
10 ^ hypothetical unobserved records N. 



the jump. In (d) we see that if there is no jump, the 
fldelity with px is always greater for ps- 

We confirm that smoothing enables better state esti¬ 
mation on average by calculating the av erag e purity, for 
10^ realisations of Y. Recall from Eq. (15) that higher 
purity means higher fldelity with the true state. We plot 
this in Eig. for two different local oscillator phases: 
4> = 7 r /2 (T-homodyne) in (a) and 4> = 0 (X-homodyne) 
in (b). Because the driving of the atom causes ay to oscil¬ 
late at a frequencies ^4 ^ 7 , it is harder to track the state 
of the atom using Y -homodyne detection, and the purity 
of the Altered state is lower than for X-homodyne detec¬ 
tion [50]. It is the former case which shows the greatest 
improvement in purity under smoothing: about 26% of 
the purity lost, because of the unobserved radiation, is 
recovered in (a) compared to about 12 % in (b). 

One can easily show m that our theory of quantum 
state smoothing includes as a special case the HMM that 
applies [inmiiis] to quantum systems that (unlike our 
atomic example) have no coherences and so are effectively 
classical. Eor genuinely quantum systems there are many 
questions about state-smoothing to explore, including: 
what happens if one assumes the unobserved unravelling 
to be different from photon counting; is there a relation 
between quantum state smoothing and the “most prob¬ 
able path” formalism of Refs. (STJ [52]; does the HMM 
inevitably emerge in the classical limit, and does quan¬ 
tum smoothing necessarily work best in that limit; and 
what experiments would show uniquely quantum aspects 
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FIG. 4: Average purity for the case considered for filtered and 
smoothed states. The average has been calculated with 10^ 
records, each one of them calculated with 10^ estimated 
^ records. The figure shows the results obtained form a 
Y-homodyne measure (top) or X-homodyne (bottom). 


of quantum state smoothing? 
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FILTERING AND QUANTUM TRAJECTORIES 
Master Equation (Unconditioned states) 


The unconditioned state’s evolution of the quantum systems considered in the paper can be described by a quantum 
master equation. There are many forms of this equation, of which the most general is the Lindblad form [1]: 

N 

p = Cp=-i[H,p\^^V[ak\p, ( 1 ) 

k=l 


where {dk} are the Lindblad operators, arising from the coupling of the reservoir to the system and whose action on 
an arbitrary state p is given by 


V[ak] • = ak 


• a 


I - 


( 2 ) 


and H = W. For example, in the particular case of two environments (a and h) interacting with the system, the 
unconditioned state’s evolution equation of the system is 


p = —i[H,p] pU cpc^ 


1 

2 


{Pb + c^c, p}. 


( 3 ) 


Quantum trajectories 

In estimation theory we are looking to describe the system conditioned upon measurements. The quantum trajectory 
approach has been shown to be analogous to the classical filtering method. In this case there is not one solution of 
the master equation, but rather different quantum trajeetories or paths, depending on the results obtained. That is, 
a quantum trajectory describes the path taken by the state of the system in the correspondent space of the density 
matrix, over time, conditioned on the measurements results Vf in each infinitesimal interval [t + dt). This process is 
described by a set of measurement operations (completely positive maps) that evolve the (unnormalized) state 

forward in time: 


- Vfr* P^^ (t) ■ (4) 

Here the average of these unnormalized states yields the master equation when an appropriate ostensible distribution 
for the result is used: 

post(rt) Mrt = 1 + Cdt, ( 5 ) 

rt 

where we consider to be discrete for simplicity. Clearly cannot be defined independently of post(i*t) or Iri 
the case of complete observation, this preserves the purity of the state, and we can write 

A4r,. = Mr,.My (6) 

Starting with p(to) = po? fhis procedure generates the state conditioned on the whole past record: (t). This is 

an unnormalized state (as indicated by the tilde) analogous to the unnormalized classical state from Eq. (2) in the 
paper. However, this state can be renormalized, to define the filtered state, with 


PF{t) = p^p) 
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The unnormalized state carries information about the actual probability distribution of the results measurements as 
follows 

p(fet|po) = Tr[p^^(t)] Post(fetlpo)- (7) 

The quantum trajectories can be discontinuous (quantum jumps) or continuous (quantum diffusive trajectories) 
depending on the kind of measurements that are conditioning the evolution. 


Diffusive Quantum Trajectories 

We consider now the case where the measurements are diffusive in nature and the record in each interval is a 
continuous variable. For example, for the case of unit-efficiency homodyne detection with a single channel 6, the 
result yt is a real number and we can choose to be a zero-mean Gaussian of variance 1/dt. Then the 

measurement operation is defined by = My^ • where [2] 

My,=l- iHdt - (8) 

is the measurement operator acting only in the system Hilbert space. It is easy to verify that this reproduces the 
master equation evolution to order dt: 

j dyt post{yt) = 1 - i[H, •] + V[h] • . (9) 

Here is the phase of the local oscillator. We say we make an X-homodyne measurement when = 0 and Y- 
homodyne when = 7r/2. 

A normalised version of the filtering state can be expanded using Ito calculus, yielding 

+V[b]p*y^{t)^ dt+ (e*‘*’(S- (S))p^^(t) + H.C.) dw{t). (10) 

Since this is a normalized state, it averages to the master equation solution according to the actual probability 
distribution, for which dw{t) = ytdt — Tr[e*^6 -h e~'^^b^]dt is a Wiener increment satisfying Ey^[dw{t)] = 0 and 
Ey^[dw‘^{t)] = dt. In this form it is elementary to see that averaging over the measurement results yields the master 
equation again. We can use this equation to directly calculate the normalised pp = P\f • 


Discrete Quantum Trajectories (Quantum Jumps) 

A similar definition, or the case of a single-channel with unit-efficiency photon counting: 

Adrift = Mrit • ^rit’ 

In this case the record is discrete — there are only two possible results in the interval [t, t dt): one photon or none. 
Choosing ostensible probabilities 

Post(^t := 1) = Adt ; post(^t := 0) = 1 - Adt, 
the correspondent measurement operators are [2] 

Ml =c/VX, 

Mo = l-(iH+ I (etc - A)) dt. 

Again, it can easily be verified that averaging by using these ostensible probabilities gives 

J2nt Post(nt) pntit) = (1 - Adi).Mo p{t)M^ + XdtMi p{t)Ml 
= p{t) — i[H^ p{t)]dt -f 'D[c]p{t)dt 

as required for consistency with the master equation. The actual probabilities, which relate to the normalized state, 
are determined by 

p{nt := 1) = postirit := l)Tr[Mip{t)Ml] = Tr[p{t)c^c]dt, 
and similarly for p{nt := 0), which evaluates to 1 — p(nt = 1). 


( 11 ) 

( 12 ) 

(13) 


( 14 ) 
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Multiple channel monitoring 

When the system is monitored through more than one channel, as we must consider for the case of the quantum 
smoothing state method, the conditioned state is defined in terms of the measurement operators of both channels. 


— ^yt ^rit • ^rit ^yt' 

This operation gives Eq. (3) when averaging over the results for both monitorings if we take the Hamiltonian in each 
measurement operator to be ^/2 rather than the full H. In the case of average over diffusive trajectories, we have to 
order dt. 


J Post{yt)-Mnt,yt* — J postivt) ^yt {^rit • ^nt)^yt 

= Mn, • Ml - l[H, •Jdt + V[b] • dt 

= Mnt • . 

and similarly for discrete measurements, 

E Post(nt) = E Post(nt) My, Mn, • 

nt fit 

= My, •Ml,-i [H, •]dt + V[c] • dt 


(16) 


(17) 


Just as for the individual jump and diffusion case, we can easily normalize the state, and generate records with the 
correct actual distribution. 


Choosing the Ostensible probabilities 


The ostensible probability distribution for the measurement results arbitrary, as long as it is fixed for a given 
ensemble for simulation purposes. Thus, for the purpose of smoothing, where we have to simulate many different 
unobserved records ^ for a given observed record O , the ostensible distribution for ^ could depend upon O (since 
this is fixed). For the most efficient simulations (that is, giving the smallest statistical errors for a given ensemble size), 
we would like the ostensible distribution p°®^(U) to be as close as possible to the actual distribution | O ). This 
last is, however, exactly what we are trying to determine from the simulations — the smoothed distribution ps(^) — 
so we cannot expect to just guess it. Moreover, the more complicated is, the more complicated the quantum 

trajectory equations corresponding to this ostensible distribution will be. Thus, for simplicity, we choose an ostensible 
probability distribution using only the past information giving a distribution that factorizes at different times. 
For the example in the paper, where O = Y and U = N, the distribution for rit is thus described by a parameter A 
as introduced in Eq. (11), but this becomes time-dependent and determined by the filtered state Pf(^) = (^)- 


X(t) = Tr[c‘*‘cpF(t)]. 


(18) 


That is, we choose the ostensible rate of unobserved jumps (photon counts) to equal the rate that would pertain, 
as in Eq. (14), if Pf(^) were the state of the system (which it is not — the unknown ^ (t) is the True’ state 
PT{t)). A more sophisticated approach, which may give somewhat smaller standard deviation, would be to use a small 
ensemble to generate an approximation to ps(t), and then use this, in place of pF(t), in Eq. (18), to generate a new 
approximation to ps(t) from a larger ensemble of unknown records. 

Note that the use of ostensible probability distributions is not limited to quantum systems. For example, the 
classical analogue of Eq. (7) is 


p(k(|po) = J dx'p^^(x() post(ktlpo) 


( 19 ) 
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SMOOTHING PROBABILITY 


In this section we will show how to manipulate the expressions relating to probabilities to find the smoothed 
probability ps(tjt) = = Prhf- = tJ^I^ ,p(to)]. Starting from the probability definition 

ps(tji)= p(tjtiht,hi) 

^ p(tjtiht,ht)p(ht,ht) p(h,ht) ^ 
p(ht, ht) p(0t,0t) 

we can use Baye’s rule to rewrite the expression as: 

p(htih,hi)p(tji,ht) 


Ps 


{ht) = 




so that, as p(A, B) = p{A\B)p{B) and p{A) = J2 b PiM^) p{B), we finally get, 

rtr 1 — P(®t|ht, Ot) p(tjt|Ot) 


( 20 ) 


AVERAGE PURITY AND AVERAGE FIDELITY 


The smoothed quantum state’s purity, and its fidelity with the true state, are appropriate measures of the model’s 
success. We will demonstrate in this section that for the case where the ‘true’ state is pure, the values of these 
two quantities, averaged over all records (observed and unobserved), is identical. This is true both for filtering and 
smoothing. 

In Eq. (9) in the paper, the quantum smoothed state is defined as an ensemble average over realizations of the 
unknown record tj^. If we define a general conditioning C, filtered and smoothed states can be described similarly, 

Pc(t) = (21) 

ti 

^^ _ 

where = O for smoothing (C=S), and = Of for filtering (C=F). From the definition for purity in Eq.(13) 


P[pc{t)]= Tr [pc{t) pc{t)] 
= Tr 


s pc(ht)p^^ (t) pc{t) 

Y, pc(tj() pc{t)]. 

V 

Also, when the true state px is pure, 

F[pc{t),pT] = Tr pc{t)pT{t)] 

= Tr Pcit) ■ 

Averaging Eq. (22) over the appropriate ensemble of measurement records O^, 

Eon{P[Pc{t)]}= E p(Oo)Ep(E|Oo)Tr[p^^_^^(t)pc(i)] 

tit 

= E) E p(E|On) p(On) Tr[/3^^ |j^(t)pc(^)] 

On tJj 

= E p(On! Ot) Tr |j^(t)pc(i) 

On,tit 

= fjd'^[^T(i),Pc(i)]}, 

where here we can think of the ensemble average as being over all possible ‘true’ records tt as well as 


( 22 ) 


(23) 


( 24 ) 
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RELATION TO CLASSICAL HIDDEN MARKOV MODELS 

In the paper we mention that the results relating to HMMs in [3, 4], and also in the prior work of [5] (which is 
actually somewhat more general), can be easily obtained as a restriction of the quantum state smoothing model. In 
these systems, the system dynamics—the quantum master equation—can be reformulated, using a suitable basis, as a 
classical master equation. That is, the stochastic dynamics can be described as transitions (quantum jumps) between 
orthogonal states of the system. This means that, independently of whatever record O an observer Alice has, a 
second, hypothetical observer (Bob) can always know exactly what state the system is in. That is, Bob’s unobserved 
(by Alice) jump record completely determines the state of the system at time t, as the rank-one projector 
The label x, identifying which state the system is in at time t, is the HMM for the system. 

We will now show in detail the connection between our quantum state smoothing technique and the HMM technique. 
Recall our method for calculating the quantum smoothed state: 

ps{t) = ^ 

where ps(^) = p(^t|^)- Now in this case Introducing the completeness relation Sx'^x = 1, 

for any value of the classical variable x, and also using 

= y]ps(a:^)na;'. 

x' X 

That is, the problem reduces to the classical smoothing problem, with the simple solution (requiring no ensemble 
averaging) 

ps{x) oc p(h|a;) p(a;|ht), 

as used by Armen et al [5] and Gammelmark et al [3, 4]. 


(26) 

(27) 


we can rewrite the smoothed state as 

Psit) = 
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